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Abstract 

This paper provides detailed comparisons of chemical reaction mechanisms of H2 applicable at high preheat 
temperatures and pressures relevant to gas turbine and particularly Alstom's reheat gas turbine conditions. 
It is shown that the available reaction mechanisms exhibit large differences in several important elementary 
reaction coefficients. The reaction mechanisms are assessed by comparing ignition delay and laminar flame 
speed results obtained from CHEMKIN with available data, however, the amount of data at these conditions 
is scarce and a recommended candidate among the mechanisms can presently not be selected. Generally, 
the results with the GRI-Mech and Leeds mechanisms deviate from the Davis, Li, Conaire, Konnov and 
San Diego mechanisms, but there are also significant deviations between the latter five mechanisms that 
altogether are better adapted to hydrogen. The differences in ignition delay times between the dedicated 
hydrogen mechanisms (0 Conaire, Li and Konnov) range from approximately a maximum factor of 2 for 
the HVair cases, to more than a factor 5 for the H2/O2/AR cases. The application of the computed ignition 
delay time to reheat burner development is briefly discussed. 
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1. Introduction 

Among the technologies for fossil fuel power pro- 
duction with CO2 capture, "pre-combustion" re- 
moval of CO2 is one of the strong candidates. Pre- 
combustion technology for stationary power pro- 
duction implies burning hydrogen or hydrogen rich 
mixtures in a gas turbine. Utilizing hydrogen in 
non-premixed burners without dilution (steam or 
nitrogen) causes unacceptable levels of NO^ emis- 
sions due to the high flame temperature To 
avoid the higher cost and maintenance that comes 
with dilution, there is an increasing interest in 
developing non-diluted premixed or partially pre- 
mixed hydrogen burners. 

Hydrogen is a particular molecule with high dif- 
fusivity and reactivity and there are challenges to 
safe and stable operation of burners utilizing non- 
diluted premixed hydrogen. One combustion tech- 
nology where premix combustion of hydrogen has 
a potential is the reheat or sequential gas turbine 
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technology. Other related applications are ramjet, 
scramjet and afterburner technologies. The reheat 
gas turbine technology was first commercialized in 
1948 by Brown Boveri Co., and is marketed today 
by Alstom. These turbines had the sequential com- 
bustor design, which is used to increase efficiency 
and provide operational flexibility while manag- 
ing low emissions. The combustion system uses 
an Environmental (EV) burner in the first com- 
bustion stage followed by a Sequential EV (SEV) 
burner, illustrated in Figure [TJ in the second stage 
[H, Q • Since the reheat combustor is fed by high- 
temperature expanded exhaust gas of the first com- 
bustor, the operating conditions allow auto-ignition 
(spontaneous ignition) of the fuel air mixture with- 
out additional energy being supplied to the mix- 
ture. To prevent ignition of the fuel air mixture 
in the mixing region, an appropriate distribution 
of the fuel across the burner exit area must be ob- 
tained and the residence time in the mixing region 
must not exceed the auto-ignition delay time. As 
of today with conventional fuels, this is solved by 
using delta wing shaped vortex generators to mix 
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Figure 1: Operation principle of Sequential Environmental 
(SEV) combustor. Illustration from Ciani et al. Q] 



the fueL which is injected in a centrally positioned 
lance 

The chemical mechanism is an important element 
in several tools for gas turbine calculations. Com- 
putational Fluid Dynamics (CFD) has evolved to 
become one of the main tools in design of gas tur- 
bine combustors, especially in the developing stage 
where experiments are relatively resource intensive 
and expensive. Accurate CFD predictions of hydro- 
gen combustion rely on high precision in the chem- 
ical mechanisms. 

A number of detailed mechanisms for hydrogen, 
often known as the core of any detailed hydrocar- 
bon combustion reaction mechanism have been 
developed during the last two to three decades. 
Some of these are sub-mechanisms of larger hydro- 
carbon mechanisms such as methane or ethane. A 
brief overview of the mechanisms considered in the 
present work is given in the following and summa- 
rized in Table [T] The H2 / O2 mechanism of Li et 
al. @ is based on the mechanism that originally 
was developed by F.L. Dryer at Princeton Univer- 
sity and further developed by Mueller et al. 0- 
This mechanism has been compared against a wide 
range of experimental conditions with temperatures 
between 298-3000 K, pressures from 0.3-87 atm and 
equivalence ratios from 0.25-5. The mechanism of 
Mueller was also used as the basis for the mecha- 
nism developed by O Conarie et al. [8j at University 
of Ireland in Galway. The series of experiments 
numerically investigated ranged from 298-2700 K, 
pressures between 0.05-87 atm and equivalence ra- 
tios from 0.2-6. Konnov [5j has derived a H2/O2 
mechanism from a methane mechanism with a re- 
cent review and update of the elementary reactions. 
This mechanism is validated with ignition experi- 
ments (950 to 2700 K, sub-atmospheric to 87 atm) 
and flame speeds at pressures from 0.35-4 atm. 
Davis et al. |9| developed a H2/CO/O2 mechanism 



based on an extensive validation against ignition 
delay and flame speed data in the range from 0.05- 
64 atm and 298-1754K (for hydrogen). The San 
Diego mechanism developed by Forman Williams 
and co-workers at the University of California, San 
Diego 1(J 11 1 is also studied here. The full mech- 
anism includes hydrocarbons up to C3 [12]. The 
H2/O2 part has been validated against hydrogen 
ignition delay data up to 33 bar [13j . These five 
mechanisms have been selected for consideration 
based on thoroughness in validation against hydro- 
gen ignition and flame speed data. In addition 
we consider two of the most widely used hydro- 
carbon mechanism that includes a H2/O2 subset. 
The Leeds methane mechanism has been developed 
at Leeds University by Hughes et al. [14|, while 
the GRI mechanism origins from the Gas Research 
Institute [15j. More comprehensive summaries of 
the various mechanisms for hydrogen combustion 
are for example given by del Alamo et al. [l3[ and 
Strohle and Myhrvold [161] . 



Mechanism 


N s 


N R 


K 


Reference 


Li 


10 


23 


23 


m 


O Conaire 


10 


21 


21 


[8]/[17] 


Konnov 


10 


33 


33 


M 


Davis 


14 


38 


25 


WUSH 


San Diego 


46 


235 


21 


[10j/[12] 


Leeds 


37 


175 


23 


M/\m 


GRI-Mech 


53 


325 


28 


[15] 



Table 1: Overview of the chemical kinetic mechanisms and 
their most recent references considered in the present work. 
N s is the total number of species in the mechanism (includ- 
ing Argon). Nr is the total number of reactions in the mech- 
anism while N R is the number of reactions in the H2/O2 
subset (duplicate reactions and separately formulated third 
body reactions are counted). 



The number of available mechanisms and vali- 
dation studies are significant, and still there are 
uncertainties related to the elementary reactions 
involved Very recently, Burke et al. [2(| in- 
vestigated the mass burning rate (the product of 
laminar flame speed and unburnt gas density) and 
found that there is a negative pressure dependence 
at high pressures and low temperatures. According 
to their results, this dependence is not captured by 
the available mechanisms due to uncertain and even 
missing third-body reactions. 

The aim of the present study is to investigate the 
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performance of the selected detailed chemical mech- 
anisms at the high temperature and pressure condi- 
tions relevant to gas turbine and particularly reheat 
engine conditions. As a starting point in reheat 
burner design, the burner residence time is chosen 
by considering the ignition times. This is based 
on the assumption that the flame position is con- 
trolled by pure auto-ignition. The effect of pressure 
on ignition delay time and laminar flame speed is 
investigated in order to illustrate the effect of vari- 
able engine load. For the universality of this study, 
pressures also higher than the baseload operating 
pressure of the Alstom SEV burner are considered. 
The data of Herzler and Naumann [2l| and Petersen 
et al. [22j are applied for comparison. There is gen- 
erally little data available for auto-ignition and lam- 
inar flame speed in the pressure range from 15 to 
30 bar particularly at high temperatures. Conse- 
quently the mechanisms are also compared at con- 
ditions where they have not been validated against 
experiments. Nevertheless, for predictions of igni- 
tion delay and flame speed under reheat gas turbine 
conditions, their performance at these conditions is 
of interest. 

2. Overview of important elementary reac- 
tions 

In Table [5] the most important H2/O2 reactions 
for ignition, extinction and flame propagation iden- 
tified by Strohle and Myhrvold [23] are listed. We 
have also verified that for the results in the present 
work this 11-step reduced Li mechanism closely re- 
produces the results with the full Li mechanism. 
This has, however, not been verified for the other 
selected mechanisms. The reactions coefficients 
marked with (*) are for the reverse reaction and can 
not be directly compared. Among the selected reac- 
tions, the reverse reaction (10) is the main initiating 
step which creates H-radicals. The chain branch- 
ing/propagating steps (1-3) are important for the 
auto-ignition, but competes with the chain termi- 
nating steps (9a-f) that suppress the formation of 
the highly reactive H-radical 0, El. The chain 
terminating steps are pressure dependent and be- 
come increasingly important at higher pressures. 

In the following, a brief overview of the differ- 
ences between the reaction coefficients listed in Ta- 
ble [2] and [3] are given. More comprehensive and 
in-dept discussions of these coefficients are found 
elsewhere in the literature [5J, [2J, |25j . The chain- 
branching reaction (1) is recognized as one of the 



most sensitive and important reactions to hydro- 
gen combustion, and the differences between the 
reaction coefficients are also significant. For ref- 
erence, Konnov [f| applies the latest expression 
recommended by Baulch et al. 25]. The coeffi- 
cients in the reaction (2-3) show little deviation 
for all mechanisms. Li et al. 6] modified the pre- 
exponential factor of reaction (8a) to improve flame 
speed predictions. Apart from this, the mecha- 
nisms deviate mainly in the third-body collision 
efficiencies. It should for instance be noted that 
the third-body effect of H2O in Davis is about half 
the quantity in the Li, O Conarie and San Diego 
mechanisms. Konnov applies a separate equation 
(8b) listed in Table [3] when the collision partner 
is H2O. This is effectively to take temperature- 
dependent third-body effects into account [13]. For 
the chain terminating reaction (9a) , the GRI-Mech 
and Leeds mechanisms do not include low-pressure 
limits. In addition, GRI-Mech and Leeds have sep- 
arate expressions for the collision partners H2/H2O 
and H2O, respectively. Also Konnov applies sepa- 
rate fall-off expressions for individual collision part- 
ners. Among the other mechanisms, the differ- 
ences in the low-pressure limit coefficients and in 
collision efficiencies especially for H2 and H2O are 
the most apparent differences. The Li-mechanism 
has an alternative reaction (9a) when the main 
bath is AR where the low pressure limit coefficients 
are [9.0T0 19 , —1.5, 492], the broadening factor is 
F c =0.5, and the third body efficiencies are 02=1.1, 
H20=16 and H2=3. The reverse reaction (10) is the 
most important initiating step in the mechanisms 
[loj . The differences between Konnov, which uses 
the coefficients recommended by Baulch et al. (25j |. 
and Li, O Conarie and San Diego should be noted. 
Konnov [5| reports large uncertainty of the HO2 
reaction (13). The difference between the mecha- 
nisms is evident in that additional duplicate reac- 
tions are introduced in the GRI-Mech, Davis and 
Konnov mechanisms. However, the coefficients in 
the respective other mechanisms are close to equal 
for this reaction. For reference, Li uses the expres- 
sion recommended by Baulch et al. 24|,|25| for reac- 
tion (14). All mechanisms except Davis use dupli- 
cate reactions with similar coefficients for reaction 
(14) and (14dup). Considering reaction (17) the 
mechanisms exhibit large variations. Again, Li, O 
Conarie and San Diego have the closest agreement, 
while Konnov, which adopts the Baulch et al. [25j 
recommendations, has lower pre-exponential factor 
and activation energy. The expressions applied in 
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GRI-Mech and Davis involve temperature depen- 
dences of the pre-exponential factor. 



Reaction 


GRI- 
Mech 


Konnov 


Leeds 


(8b) 

H+OH+H 2 
= H2O+H2O 




1.0-10 26 « 

-2.4 

120160 




(9b) H+0 2 (+0 2 ) 
= H0 2 (+0 2 ) 




4.7-10 12 
0.4 0.0 

[5.7T0 18 
-1.1 0.0] 
F c =0.5 




(9c) 

H+0 2 (+H 2 0) = 
H0 2 (+H 2 0) 




9.110 12 
0.2 0.0 
[3.7-10 19 
-1.0 0.0] 
F c =0.8 




(9c) 

H+0 2 (+AR) 
= H0 2 (+AR) 




4.610 12 
0.4 0.0 
[7.4-10 18 
-1.2 0.0] 
F c =0.5 




(9e) H+0 2 +H 2 
= H0 2 +H 2 


1.1310 19 
-0.8 0.0 




6.8910 15 
0.0 -2076 


(9f) H+0 2 +N 2 = 
H0 2 +N 2 


2.6-10 19 
-1.2 0.0 






(15b) 

OH+OH(+H 2 0) 
= H 2 2 (+H 2 0) 




1.010 14 
-0.37 0.0 
[1.45-10 18 
0.0 0.0] 
F c =0.5 





Table 3: Overview and comparison of reaction rates specific 
to the GRI-Mech, Konnov and Leeds mechanisms. See the 
caption of Table [2] for explanation. 



3. Method 

Ignition delay time is computed by solving the 
perfectly stirred closed reactor equations 



dt 



dh 
~dt 



0. 



dp 
~dt 



0, 



(1) 



where Yi and w,; are mass fraction and chemical 
reaction rate of species i, respectively. The initial 
value problem is solved at constant enthalpy h (adi- 
abatic) and pressure p conditions using an in-house 
code where the CHEMKIN library 26] is used for 
the property and source term calculations. Un- 
less otherwise stated, the ignition delay time is de- 
fined by the time when the temperature has reached 
200 K above the initial value. A similar criterion is 
used by e.g. Konnov 

The freely propagating one-dimensional laminar 
premixed flame calculations were performed using 



the PREMIX code [27|. The effect of species trans- 
port by temperature gradients, the Soret effect, and 
multicomponent diffusion is included in the compu- 
tations since these phenomena are of particular im- 
portance in hydrogen combustion. The first-order 
upwind discretization scheme was used with suf- 
ficient grid refinement so that a grid-independent 
solution was achieved for all cases. Since the pre- 
mixed reactant mixture is above its auto-ignition 
temperature it is important to ensure that the flame 
is positioned close enough to the inlet boundary so 
that the results are not affected by auto-ignition. 
On the other hand the flame must not be positioned 
too close, since the flame speed calculations are as- 
suming zero gradient on the inlet boundary. Due 
to these high-temperature related restrictions, the 
flame speed results were not completely indepen- 
dent of the model input fixed temperature Tfix0 
A variation in Tpix with ±50 K, which generally 
was chosen close to the inlet temperature value, 
resulted in a variation in flame speed with 4-5 %. 
Under these high temperature and pressure con- 
ditions, achieving experimental results of laminar 
flame speeds will be extremely challenging. Despite 
of this the laminar flame speeds have been included 
in this work because generally in burner design both 
ignition delay and flame speed are used. 

4. Results and Discussion 

4--1. Ignition delay results 

Petersen et al. [22| measured the ignition de- 
lay times of stoichiometric H2/O2 mixtures using 
a high-pressure shock tube [l^ . The measurement 
data considered in the present work were performed 
at 33 bar between approximately 1175 K and 1300 K 
with 2% H 2 , 1% 2 and 97% dilution with Ar- 
gon. Due to the relatively low temperature in- 
crease, the criteria for determining the ignition time 
in the computations was set to when the temper- 
ature reached 100 K above the inlet temperature. 
Figure [5] show that the GRI-Mech and the Leeds 
mechanism largely overpredict the ignition delay 
data of Petersen. The latest Li mechanism in- 
volves a compromise on the equation (9a) where 
the low pressure limit and third body efficiencies 
arc different depending on whether the main bath 



x When solving the freely propagating flame problem, an 
additional boundary condition is required. This condition is 
the temperature Tpix of the flame at a given position 27]. 
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Reaction 


GRI- 

Mech 


Davis 


Konnov 


Li 


O Conaire 


San Diego 


Leeds 


(1) H+O2 = OH+O 


2.6510 16 

-0.7 

17041 


2.6410 16 

-0.7 

17041 


2.06-10 14 

-0.1 

15022 


3.55-10 15 

-0.4 

16599 


1.9110 14 
0.0 16440 


3.5210 16 

-0.7 

17070 


9.76-10 13 
0.0 14821 


(2) 0+H 2 = OH+H 


3.87-10 4 
2.7 6260 


4.59-10 4 
2.7 6260 


5.06-10 4 
2.7 6290 


5.08-10 4 
2.7 6290 


5.08-10 4 
2.7 6292 


5.06-10 4 
2.7 6291 


5.12-10 4 
2.7 6277 


(3) Ha+OH 
H2O+H 


2.1610 s 
1.5 3430 


1.7310 s 
1.5 3430 


2.14-10 8 
1.5 3450 


2.16-10 8 
1.5 3430 


2.1610 s 
1.5 3430 


1.1710 9 
1.3 3635 


4.52-10 8 M 
1.6 18401 


l^oaj xl-f-L-'rli-lVl — 

H 2 o+M 


on 1 n22 

-2.0 0.0 

H 2 0=3.65 

H 2 =0.73 

AR=0.38 


A An 1 f)22 
Q.QV) - 1U 

-2.0 0.0 

H 2 0=6.3 

H 2 =2.0 

AR=0.38 


« nfi 1 n27 (*) 

-3.3 

120770 
H 2 O=0 
H 2 =3.0 
N 2 =2.0 
2 =1.5 


q on t n22 

-2.0 0.0 

H 2 0=12 

H 2 =2.5 

AR=0.38 


A crn 1 n22 

-2.0 0.0 

H 2 0=12 
H 2 =0.73 
AR=0.38 


a nn 1 n22 

4.UU- 1U 

-2.0 0.0 

H 2 0=12 

H 2 =2.5 

AR=0.38 


^ t n22 

-2.0 0.0 

H 2 0=2.55 
N 2 =0.4 
O 2 =0.4 
AR=0.15 


(9a) H+0 2 (+M) = 
H0 2 (+M) 


2.8010 18 
-0.9 0.0 
U2— u 
H 2 0=0 
N 2 =0 
AR=0 


5.1210 12 

0.4 0.0 

[O.o-lU 

-1.4 0.0] 

F c =0.5 

O 2 =0.85 

H 2 0=12 

H 2 =0.75 

AR=0.4 


4.66-10 12 

0.4 0.0 

re 7 1 r>19 
[0. 1 *1U 

-1.4 0.0] 

F c =0.5 

O 2 =0 

H 2 0=0 

H 2 =1.5 

AR=0 


1.48-10 12 

0.6 0.0 

["fi A 1 n20 
[D.4- 1U 

-0.17 

5248] 

F c =0.8 

O 2 =0.78 

H 2 0=11 

H 2 =2 


1.48-10 12 

0.6 0.0 

[o.o-lU 

-0.4 

-1120] 

F c =0.5 

O 2 =0.78 

H 2 0=14 

H 2 =1.3 


4.6510 12 
0.4 0.0 

[0. I'W 

-1.4 0.0] 

F c =0.5 

H 2 0=16 

H 2 =2.5 

AR=0.7 


2.10-10 18 
-0.8 0.0 

f \ r» A 

\J2 — U.4 

H 2 0=0 
H 2 = .75 
N 2 =0.67 
AR=0.29 


(10) HO2+H 
H2+O2 


4.48-10 13 
0.0 1068 


5.92-10 5 M 
2.4 53502 


1.05-10 14 
0.0 2047 


1.66-10 13 
0.0 823 


1.6610 13 
0.0 820 


1.6610 13 
0.0 823 


4.28-10 13 
0.0 1408 


(11) HO2+H 
OH+OH 


8.40-10 13 
0.0 635.0 


7.48-10 13 
0.0 295.0 


1.90-10 14 
0.0 875.0 


7.08-10 13 
0.0 295.0 


7.0810 13 
0.0 300.0 


7.0810 13 
0.0 295.0 


1.69-10 14 
0.0 883.0 


(13) HO2+OH = 
H2O+O2 


1.45-10 13 
0.0 -500 


2.3810 13 
0.0 -500 


2.89-10 13 
0.0 -500 


2.89-10 13 
0.0 -497 


2.8910 13 
0.0 -500 


2.8910 13 
0.0 -497 


2.89-10 13 
0.0 -500 


(13dup) HO2+OH = 


5.0010 15 

U.U 1 1 OOU 


1.0010 16 

U.U 1 1 OOU 


9.27-10 15 

U.U 1 / ouu 










(14) HO2+HO2 = 
H2O2+O2 


4.20-10 14 
0.0 12000 


3.6610 14 
0.0 12000 


1.03-10 14 
0.0 11040 


4.20-10 14 
0.0 11982 


4.2010 14 
0.0 11980 


3.0210 12 
0.0 1386.2 


4.22-10 14 
0.0 11957 


(14dup) HO2+HO2 = 
H2O2+O2 


1.3010 11 
0.0 -1630 


1.3010 11 
0.0 -1630 


1.94-10 11 
0.0 -1409 


1.30-10 11 
0.0 -1629 


1.3010 11 
0.0 -1629 




1.32-10 11 
0.0 -1623 


(15a) OH+OH(+M) 
= H 2 2 (+M) 


7.4010 13 

-0.4 0.0 

[2.3-10 18 

-0.9 

-1700] 

H 2 0=6 

H 2 =2 

AR=0.7 


1.1110 14 

-0.4 0.0 

[2.0-10 17 

-0.58 

-2293] 

H 2 0=6 

H 2 =2 

AR=0.7 


1.00-10 14 
-0.4 0.0 
[2.4-10 19 
-0.8 0.0] 
F c =0.5 
H 2 O=0 


2.95-10 14 (*> 
0.0 48430 
[1.2-10 17 
0.0 45500] 
F c =0.5 
H 2 0=12 
H 2 =2.5 
AR=0.64 


2.95-10 14 <*) 
0.0 48400 
1.3-10 17 
0.0 45500] 
F c =0.5 
H 2 0=12 
H 2 =2.5 
AR=0.64 


7.4010 13 

-0.4 0.0 

[2.3-10 18 

-0.9 

-1700] 

H 2 0=6 

H 2 =2 

AR=0.4 


7.23-10 13 
-0.4 0.0 
[5.6-10 19 
-0.76 0.0] 
H 2 0=6.5 
N 2 =0.4 
O 2 =0.4 
AR=0.35 


(17) H2O2+H = 
H0 2 +H 2 


1.21-10 7 
2.0 5200 


6.05-10 6 
2.0 5200 


1.70-10 12 
0.0 3755 


4.82-10 13 
0.0 7950 


6.03- 10 13 
0.0 7950 


4.79-10 13 
0.0 7959 


1.69-10 12 
0.0 3747 



Table 2: Overview and comparison of reaction rate coefficients for the most important H2/O2 reactions. For each mechanism 
and reaction the pre-factor, temperature exponent and activation energy are listed in respective order. The low-pressure limits 
of the reactions arc given within brackets [ ]. The unit of activation energy is cal/mole. Numbering of the reactions is in 
accordance with Li et al. [f|. Reactions marked with (*) are reverse reaction coefficients. 
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is N2/He or Ar. This is done in order to take 
into account differences in broadening factors and 
temperature-dependences of collision efficiencies for 
different bath gases [6(. When the "wrong" option 
is used, the Li mechanism end up on approximately 
the GRI-Mech and Leeds level. The Li, O Conairc 
and San Diego mechanisms give satisfactory agree- 
ment with the experiments in the lower tempera- 
ture range and these three mechanisms also follow 
each other very well. At higher temperatures Kon- 
nov and Davis provide a slightly closer match with 
the experiments than the respective other mecha- 
nisms. The overprediction of the experiments at the 
highest temperatures ranges from approximately a 
factor two for Konnov and Davis to more than a 
factor 6 for GRI-mech. 

Recently, Herzler and Naumann [2l[ investigated 
the ignition of methane/ethane/hydrogen mixtures 
with H 2 content from % to 100% at 1, 4 and 
16 bar. Only the measurement data with 100% H2 
are considered in the present work. The H2/O2 
mixture is diluted by 93.0% and 91.1% Argon on 
mole basis at 0=0.5 and 0=1, respectively. The 
data at 16 bar arc compared with ignition delay 
computations at 0=0.5 in Fig. [3] and at 0=1 in 
Fig. |U The GRI-Mech and Leeds mechanisms seem 
to give the highest overprediction of ignition delay 
at all temperatures. The best agreement at low 
temperatures are found for the Davis and Konnov 
mechanisms, while Li, O Conaire and to some ex- 
tend San Diego provide a closer match with the 
experiments at higher temperatures. 

Data from Herzler and Naumann are also selected 



Figure 3: Ignition delay time results with H2/O2 diluted by 
93.0% AR at 16 bar and <j>=0.5. Comparison with shock tube 
data of Herzler and Naumann |2 li . 




lOOO/Tp- 1 ] 

Figure 4: Ignition delay time results with H2/O2 diluted by 
91.1% AR at 16bar and <j>=l. Comparison with shock tube 
data of Herzler and Naumann l2li . See Fig. [3] for legend. 



6 



ioV 
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Figure 5: Ignition delay time results showing the pressure 
effect for the composition given in Fig.|4]at T=1060 K. Com- 
parison with shock tube data of Herzler and Naumann [2l| 
at 1, 4 and 16bar. See Fig.[3]for legend. 
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Figure 6: Ignition delay time results showing the pressure 
effect for the composition given in Fig.|4]at T=1160K. Com- 
parison with shock tube data of Herzler and Naumann |2 j| 
at 1, 4 and 16 bar. See Fig. [3]for legend. 



to illustrate and compare the pressure effect on ig- 
nition delay time. The data are selected at approx- 
imately 1060 K and 1160 K@ and compared to com- 
putations in Fig.[S]and Fig. El All mechanisms cap- 
ture the trend of decreasing followed by increasing 
ignition delay time with increasing pressure. The 
shift of the minimum ignition delay time to higher 
pressures for higher temperatures is observed for 
all mechanisms. This effect is explained for the 
H2-air cases in the following paragraph. The lim- 
ited amount of data points is insufficient in order to 
draw any conclusions on the performance of the dif- 
ferent reaction mechanisms. The Konnov and Davis 
mechanisms seem to predict the pressure effect best 
at 1060 K, while the Li and O Conaire mechanisms 
are closer to the experiments at 1160 K. It is worth 
noticing that the difference between the dedicated 
H2 mechanisms (O Conaire, Li and Konnov) is close 
to a factor 5 (around 8 bar) for the T=1060K case 
and approximately a factor 2.5 (around 16 bar) for 
T=1160K. 

The remaining investigation of ignition delay 
time is performed with H2-air at pressures and tem- 
peratures relevant to gas turbine and reheat engine 
conditions. No experimental work is available in 
the literature at these conditions so this part is lim- 
ited solely to a comparison between the individual 
mechanisms. The Figs. 1^1 and [TU1 show the ignition 



2 The data are selected at 1061 K, 1064K, and 1059K at 
1, 4, and 16 bar in Fig. [5] and at 1153 K, 1159 K and 1160 K 
at 1, 4, and 16 bar in Fig. [6] 



delay plotted against pressure for an initial H2-air 
mixture at a fuel-air ratio of $=0.2 and 1173 K and 
1273 K, respectively. The following 4 figures cover 
the same temperature range at $=0.75 and $=1.5. 
The pressure trend, which is similar to the trend in 
the Figs. [5] and |H1 is captured by all mechanisms. 
The ignition delay time drops going from 2 bar to a 
minimum which lays in the range 5-20 bar and then 
rises again in the interval 10-30 bar. A sensitivity 
study performed for the Li-mechanism at $=0.75 
and 1273 K showed that the chain branching re- 
action (1) dominates at low pressures. At higher 
pressures the main chain terminating reaction (9) 
becomes increasingly important and the ignition de- 
lay time decreases. Similar sensitivity is expected 
for the other mechanisms. With increasing temper- 
ature, the minimum inflection point of the ignition 
delay is predominantly shifted to higher pressures. 
According to e.g. Herzler and Naumann 21], this 
effect is due to the activation energy of reaction (1) 
which increases the rate of reaction (1) relative to 
reaction (9) with increasing temperatures. 

The discrepancies between the mechanisms are 
generally largest at higher pressures. This is ex- 
pected since the mechanisms are better validated at 
lower pressures. The GRI-Mech and Leeds mecha- 
nisms predict higher ignition delay times than the 
other mechanisms at all temperatures, pressures 
and fuel-air ratios considered. The overprediction 
of the ignition delay time with Leeds and GRI-Mech 
may be attributed to the differences in the chain 
terminating reaction (9a). GRI-Mech and Leeds 
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Figure 9: Ignition delay results with H2-air at <E>=0.75 and 
T=1173K. 



Figure 7: Ignition delay results with H2-air at <3>=0.2 and 
T=1173K. 
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Figure 8: Ignition delay results with H2-air at $=0.2 and 
T=1273K. See Fig. [7] for legend. 



Figure 10: Ignition delay results with H2-air at $=0.75 and 
T=1273K. See Fig.|9]for legend. 



have no distinction between the high and low pres- 
sure limits, and in addition the pre-factor is sig- 
nificantly higher. An attempt was made where re- 
action (9e) of the Li mechanism was replaced by 
reaction (9a), (9d) and (9c) of GRI-mech for the 
case in Fig. This gave an increase in ignition 
delay time approaching and almost coinciding with 
the GRI-Mech results at higher pressures. 

For $=0.2 (see Figs. [7] and © the behavior of 
the mechanisms are largely divided in three groups. 
GRI-Mech and Leeds are significantly overpredict- 
ing the main group consisting of Davis, Li, O 
Conaire and San Diego, while Konnov predicts 
somewhat shorter ignition delay times than the 
main group. When $=0.75 (see Figs. M an d ITUf 
the three groups are GRI-Mech and Leeds with the 
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Figure 11: Ignition delay results with H2-air at <I>=1.5 and 
T=1173K. 
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Figure 12: Ignition delay results with H2-air at <3?= 1 .5 and 
T=1273K. See Fig. [TT] for legend. 



longest ignition delay times, Konnov, O Conaire 
and Davis with the shortest, while Li and San 
Diego follow each other closely somewhat above the 
second group. Generally, the deviation between 
the mechanisms increases with increasing equiva- 
lence ratio, and for $=1.5 and temperatures 1273 K 
(see Figs. [Til and [T2")) . the ignition delay times are 
broadly distributed. It is interesting to note that, 
disregarding GRI-Mech and Leeds, the Davis mech- 
anism predicts the longest ignition delay times at 
lean conditions, while it predicts the fastest ignition 
at rich conditions. This may be attributed to differ- 
ences in third-body efficiencies as seen in Table [21 
Another interesting observation is that the Li and O 
Conaire mechanisms, which share the same inheri- 
tance in the Mueller mechanism 0], exhibit signif- 
icant differences except at $=0.2. The Li and San 
Diego mechanisms are very close in behavior at all 
conditions considered, even though the coefficients 
of for instance the important reactions (1) and (9a) 
are significantly different. The difference in ignition 
delay times between the dedicated hydrogen mech- 
anisms (O Conaire, Li and Konnov) ranges from 
approximately a maximum factor of 1.7 at $=0.2 
to more than a factor 2 at $=1.5. As a starting 
point in reheat burner design, the burner residence 
time is estimated from the ignition delay time. If 
the design of the reheat combustor is based on these 
computations, it implies that there is at least a fac- 
tor 2 uncertainty margin in the size of the equip- 
ment. Taking the differences of the H2/O2/AR case 
shown in Figs. [5] and [5] into account, the uncertainty 
factor may even be closer to 5. 

4-.2. Laminar flame speed results 

The laminar flame speed Sl computed with the 
PREMIX code@ are given in Figures [T2] through 
\W\ for fuel air ratios of $=0.2 and $=0.75, and 
temperatures of T=1173K and 1223 K. At these 
conditions, the H2-air mixture is above the auto- 
ignition temperature. The shortest ignition delay 
time for $=0.75 and T=1223K is approximately 
6 /is. When the flame speed at the same conditions 
is about 3800 cm/s, the distance from the upstream 
inlet at where auto-ignition occurs is 0.023 cm. 
Considering the temperature profile in the corre- 
sponding flame speed calculations (not shown), the 



3 For laminar flame calculations with the Davis mecha- 
nism an adjusted version of the CHEMKIN transport code 
must be used. The code was obtained from the developers 
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preheat-zone of the flame is located about 0.008 cm 
from the upstream inlet. This should indicate 
that the flame propagation is "ahead" of the auto- 
ignition, and the results should not be significantly 
disturbed by this. However, the highest tempera- 
ture and equivalence ratios considered in the igni- 
tion delay studies were omitted here due to prob- 
lems with convergence. We also refer to the discus- 
sion about the Tprx parameter in Sec. [3] 

Between approximately 5 atm and 30 atm the 
flame speed decreases with increasing pressure for 
all mechanisms. This is in accordance with flame 
speed theory, where the relative importance of the 
third order chain terminating reactions (9a-f) in- 
creases with increasing pressure. 

The flame speed predicted by the GRI-Mech and 
Leeds mechanisms are generally significantly lower 
than the respective other mechanisms. The excep- 
tion is the Davis mechanism which underpredicts 
the Leeds mechanism at high pressures for $=0.75. 
At T=1223K and for both $=0.2 and $=0.75 the 
flame speed slightly increases from 2 bar before it 
decreases again. This behavior is predicted by all 
mechanisms except the GRI-Mech and the Leeds 
mechanisms. A modified Li mechanism was con- 
structed where reaction (9e) was replaced by reac- 
tion (9a), (9d) and (9e) of GRI-Mech. This modi- 
fied Li mechanism did not predict the increase in 
flame speed at lower pressures, and approached 
the GRI-Mech results at higher pressures. Hence, 
some of the differences between the flame speed pre- 
dictions with the GRI-Mech/Leeds and the other 
mechanisms are attributed to the differences in the 
chain-terminating reactions (9a-f). Here the differ- 
ences between the dedicated hydrogen mechanisms 
(O Conaire, Li and Konnov) are less pronounced 
than for the ignition delay calculations. The largest 
deviation of about 30 % is found between O Conaire 
and Konnov at high pressures and $ = 0.2. 

5. Summary and conclusions 

Accurate predictions of hydrogen combustion 
rely on high precision in the chemical mechanisms. 
In the present work, the performance of selected de- 
tailed hydrogen chemical mechanisms at the high 
temperatures and pressures relevant to gas turbine, 
and particularly reheat engine conditions, are in- 
vestigated. For the universality of this study, pres- 
sures also higher than the baseload operating pres- 
sure of the Alstom SEV burner are considered. The 




Figure 13: Laminar flame speed results with H2-air at $=0.2 
and T=1173K. See e.g. Fig.[l6]for legend. 




Figure 14: Laminar flame speed results with H2-air at $=0.2 
and T=1223K. See e.g. Fig.QUfor legend. 
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Figure 15: Laminar flame speed results with H2-air at 
$=0.75 and T=1173K. See e.g. Fig. fl6lfor legend. 
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Figure 16: Laminar flame speed results with H2-air at 
$=0.75 and T=1223K. 



lated to the design of actual combustion equipment, 
as for instance the SEV burner by Alstom, the im- 
pact on cost is significant. Hence, to reduce the 
uncertainties in predictions with hydrogen mecha- 
nisms at these high temperature and pressure con- 
ditions more experiments, particularly for ignition 
delay time, are needed. 

Generally in burner design, ignition delay and 
flame speed together are needed to define the char- 
acteristics. This motivates including the laminar 
flame speed results even at the high temperatures 
considered in the present work. Future work should 
address the role of ignition time and flame speed 
towards SEV flame stabilization. Finally, basic 
knowledge of the turbulent flame speed at such con- 
ditions is crucial for the development of hydrogen 
combustion engines. 



effect of pressure on ignition delay time and lami- 
nar flame speed is investigated in order to illustrate 
the effect of variable engine load. Very little ex- 
perimental data are available at these conditions so 
a recommended candidate among the mechanisms 
can presently not be selected. The GRI-Mech and 
Leeds mechanisms seem to give the highest over- 
prediction of the ignition delay data at all temper- 
atures. The best agreement at low temperatures are 
found for the Davis and Konnov mechanisms, while 
Li, O Conaire and to some extend San Diego pro- 
vide a closer match with the experiments at higher 
temperatures. 

Generally, the results with the GRI-Mech and 
Leeds mechanisms deviate strongly from the Davis, 
Li, O Conaire, Konnov and San Diego mechanisms, 
but there are also significant deviations between the 
latter five mechanisms that altogether are better 
adapted to hydrogen. The differences in ignition 
delay times between the dedicated hydrogen mech- 
anisms (O Conaire, Li and Konnov) range from ap- 
proximately a maximum factor of 2 for the H2-air 
cases, to more than a factor 5 for the H2/O2/AR 
cases. For laminar flame speed the differences be- 
tween the dedicated hydrogen mechanisms are less 
pronounced, with the largest deviation of about 
30 % found between O Conaire and Konnov for H2- 
air at high pressure. It is shown through compari- 
son in Table [5] and [3] that the mechanisms exhibit 
large differences in several important elementary re- 
action coefficients. 

As a starting point in design, the burner residence 
time is chosen based on the ignition delay times. 
If the predicted differences in ignition time are re- 
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